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Abstract 


The mechanism of the selective dispersion of single-walled carbon nanotubes (CNTs) by 
polyfluorene polymers is studied in this paper. Using extensive molecular dynamics sim¬ 
ulations, it is demonstrated that diameter selectivity is the result of a competition between 
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bundling of CNTs and adsorption of polymers on CNT surfaces. The preference for certain 
diameters corresponds to local minima of the binding energy difference between these two 
processes. Such minima in the diameter dependence occur due to abrupt changes in the CNT’s 
coverage with polymers and their calculated positions are in quantitative agreement with pre¬ 
ferred diameters, reported experimentally. The presented approach defines a theoretical frame¬ 
work for the further understanding and improvement of dispersion/extraction processes. 

Since their diseovery, single-walled earbon nanotubes (CNTs) have attraeted major researeh in¬ 
terest due to their extraordinary meehanieal, ehemieal and eleetronie properties.^ They are metallic 
or semiconducting depending on their ehirality and as-synthesized material is normally a mixture 
of both types. For many applieations, however, purified samples of only a eertain type are in high 
demand. Purified semieondueting tubes are required, for instanee, to aehieve a large on/off ratio 
and high carrier mobility in thin-film field effeet transistors^^ and high power eonversion effi- 
eieney for photo voltaios.^^ Moreover, for optoeleetronie applieations working in a speeifie wave¬ 
length range, the sorting of semiconducting CNTs according to diameter is of great importance. 

In view of sueh demands, methods for the selective synthesis of CNTs of a eertain eleetronie 
type or ehiralites have been developedA low-eost mass produetion of seleeted CNTs is yet 
to be aehieved, however, and post-synthesis methods are often relied on.l^ A promising post¬ 
synthesis seleetion method diseovered reeently is based on the physisorption of polymers on the 
surfaee of CNTs, which has the advantage of leaving the eleetronie properties of CNT nearly un¬ 
perturbed.^^ There is a relatively long history of using polymers to disperse CNTs in aqueous or 
organic solutionsA recent finding is that, by using suitable polymers, CNTs can be selectively 
dispersed either for a specific diameter range or for certain ehiral angles.^ Among those tested, the 
;r-eonjugated polymer group of polyfluorene derivatives shows the ability to seleetively disperse 
semi-conducting CNTs.^^^^^^ In particular, the di-oetyl substituted polyfluorene (PFO) used with 
toluene as solvent prefers to disperse small-diameter semi-eonducting nanotubes with chiral angles 
bigger than about 20 degrees.!^ With longer side-chains, larger-diameter tubes ean be dispersed but 
the ehiral angle preferenee is gradually lost.l^*^ More recently, eopolymers of polyfluorene with 
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anthracene or pyridine groups were found to seleetively disperse large-diameter semieondueting 
nanotubes with high purity and high yields This fits well with the need to fabrieate photoelee- 
tronie deviees working in the infrared wavelength range.^ Large-diameter nanotubes also benefit 
from a diminishing eontaet resistanee and higher earrier mobilities. Purified semieondueting CNTs 
have been used to fabrieate high-performanee field-effeet transistors with high earrier mobilities 
and large on-off ratios.^^ 

Intensive experimental and numerieal works have been undertaken to study the eonformation 
of polymers adsorbed on the CNT surfaee. For DNAs and some bio-maeromoleeules,l^whieh pre¬ 
fer to take a helieal eonformation even in the free state, a helieally wrapped eonfiguration on CNT 
was naturally expeeted. Studies on the adsorption eonformation of linear eonjugated polymers 
are less eonelusive. For instanee, poly(aryleneethynylene)s (PAE) were found to align linearly 
along the CNT when dispersed with toluene.^ The similar poly(p-phenyleneethynylene) poly¬ 
mer, poly[p-2,5-bis(3-propoxysulfonieaeidsodiumsalt)phenylene]ethynylene, was found to form 
helieally wrapped struetures in aqueous dispersion.l^ Imaged via SEM, it was shown that when 
dispersed in ehloroform, poly(3-hexyl-thiophene) (P3HT) forms a helieally wrapped strueture on 
the surfaee of multi-walled CNTs.^ Reeently, regioregular poly(3-alkylthiophene)s (rr-P3ATs) 
were used with toluene as solvent to enrieh semi-eondueting CNTs, and MD simulations showed 
that P3ATs take quasi-linear eonformations, as adsorbed on CoMoeat CNTs.^^In 2014 Shea et al. 
reported the first experimental study on the adsorption eonfiguration of PFOs on CNTs by using 
photolumineseenee energy transfer and anisotropy measurements Their data, however, are open 
to interpretation.^ 

In the past, efforts were also made to theoretieally explain the seleetion meehanism. For DNAs, 
the intrinsie helieal nature was believed to play a erueial role in their seleetive adsorption on 
CNTs .1^ For aromatie polymers, Nish et a/.^ found that PFOs on CNT surfaees form n-fold sym- 
metrie struetures with their baekbones aligned along the tube axis. The magnitude of the binding 
energy between CNTs and polymers was shown to inerease with the tube diameter, a trend that was 
later eonfirmed by several authors If the stability of adsorption eomplexes, as indieated by 
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the binding energy, would determine the dispersibility of CNTs, the above results^^^^^^^^^l^ would 
imply that large-diameter CNTs are more easily dispersed than small-diameter ones. This, how¬ 
ever, is in eontradietion to experimental observations that PFO prefers to disperse small-diameter 
[I5li8|23| Furthermore, helieally wrapped PFO struetures on CNTs were used to explain the 
ehirality preferenee of PFO.^^ We will show below, however, that sueh helieal struetures are not 
dynamieally stable. Reeently, a eoarse-grained model was developed and used together with statis- 
tieal meehanie arguments to explain the diameter preferenee of several pyridine-eontaining eopoly- 
mers.l^But it is unelear how well the method ean be transferred to other systems. Despite these 
advanees, it is fair to say that a thorough understanding of the diameter and ehirality seleetivity of 
the polymer adsorption method is still laeking. 

Results and Discussion 

This artiele foeuses on understanding the diameter seleetivity of the polymer adsorption method 
sinee the band gap and related eleetronie/optieal properties of semieondueting CNTs are mainly 
determined by the diameter.l^ In partieular, we propose that diameter seleetivity results from a 
eompetition between the adsorption of polymers on the CNT surfaee and the bundling of indi¬ 
vidual CNTs (see Fig. [^. Our results on four relevant polymers are in exeellent agreement with 
experimentally observed diameter preferenees^ ^^ 1^1—I and, thus, resolve a eontroversy on the na¬ 
ture of the meehanism that underlies the diameter seleetion proeess. Despite the eomplexity of 
the eompetitive dispersion of CNTs, the sueeess of our simple energetie model regarding diameter 
seleetivity relies on its eorreet representation of some key faetors ineluding sterie effeets/eoverage. 

Simulations: To study the CNT-dispersion proeess, we performed elassieal moleeular dy- 
namies (MD) simulations using foree fields. Tip sonieation treatment is known to generate high, 
loeal energy densities that break bundles into individual CNTs.For the dilute polymer eon- 
eentrations used in typieal dispersion proeesses, the polymers exist as individual moleeules.^ 
Therefore, isolated, individual CNTs and polymers were assumed as the initial eonfiguration. 
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Figure 1: The proposed mechanism of diameter-selective dispersion: a competition between 
bundling of carbon nanotubes (CNTs) and adsorption of polymers on the surface of CNTs. The 
initial state of the process, given by individual CNTs and individual polymers, is a transition state 
(on the top of a potential energy hill) created by sonication. 

Solvent molecules of toluene were usually not explicitly included here. We tested that their in¬ 
clusion did not significantly change the results but mostly slowed down the adsorption dynam¬ 
ics. Four representative types of polymers were considered in this study: the homopolymer of 
polyfluorene with side-chain length C8 (PFO) or C6 (PFH), and copolymers with anthracene group 
poly[(9,9-dihexylfluorenyl-2,7-diyl)-co-(9,10-anthracene)] (PFH-A), or pyridine groups poly[9,9- 
didodecylfluorene-2,7-diyl-alt-pyridine-2,6-diyl] (PFD-Py). The chemical structures of these four 
polymers are presented in the Supporting Information. Furthermore, 13 CNTs with diameters in 
the range from 0.8 to 1.4 nm were considered. Such diameters are typically obtained in high- 
pressure CO conversion (HiPco) or pulsed laser vaporization (PLV) synthesis. Further information 
on our simulation methods can be found in the Method section. 

Adsorption complexes: The geometries of adsorption complexes were obtained by MD simu¬ 
lations using many different initial configurations, temperatures, and CNT diameters. The simula¬ 
tions always lead to an almost linear alignment of PFO chains on the CNT surface, even after using 
initial conditions that promote the formation of helically wrapped structures. Through geometry 
optimizations, we found that a multitude of such helically wrapped structures,!^ with different 
pitches and surface coverages, are local minima on the potential energy landscape (see Fig.[^(a) 
and S2).!^ However, if they were subject to MD simulations unwrapping proceeds gradually and 
after a sufficiently long run, a linearly aligned structure, as shown in Fig. |^(b), was always ob¬ 
tained. We conclude that helically wrapped adsorption complexes are metastable.!^ 
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Figure 2: The geometry of adsorption complexes: (a) A helically wrapped configuration is 
metastable, i.e., a local minimum on the potential energy landscape, (b) A snapshot of the much 
more stable, linearly aligned configuration of PFO on a (8,6) CNT. 

Binding energy and stability of adsorption complexes: A standard measure used to charac¬ 
terize the stability of the adsorption complex is the binding energy. It is defined as the difference 
between the potential energy of an adsorption complex and the sum of its constituent molecules. 
For the adsorption of polymers on CNT, it reads 

^ CNT-Polymer ^CNT -Polymer ^CNT ^Polymer- (J^l 

The binding energy for the adsorption of a single polymer chain on a CNT is shown in Fig. |^(a). 
Note that the magnitude of the binding energy increases with the tube diameter in agreement with 
previous results .^1 1 8 | 20 | 23 | caused by a better contact between polymers and CNT owing to 

the increasingly flatter surface of large-diameter CNTs.^^The side chain contributes a large part, 
about two-thirds, to the total binding energy of PFO. Consistently, the binding energy for PFH 
is smaller due to a shorter side-chain length magnitude of the binding energy of PFH- 
A is smallest, which means that, for all the tested polymers with similar length, it is the easiest 
to remove from a CNT surface. This is in qualitative agreement with our recent experimental 
observation that PFH-A can be washed away from thin films deposited using dispersed CNTs 
(unpublished results). In contrast, PFO cannot be washed away in the same manner. 

Surface coverage of CNTs by polymers and binding energy of CNT-polymer complexes: 
To avoid the rebundling of CNTs after sonication, it is necessary to sufficiently cover the CNT 
surface with polymers. 
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Figure 3: The binding energy E^^T^Poiymer adsorption complexes for (a) a single poly¬ 

mer chain and (b) the maximal coverage of the CNT surface. The magnitude of the binding energy 
increases with the nanotube diameter. Therefore, the binding energy ^^one cannot 

explain why polymers selectively disperse CNTs with specific diameters. Note the discontinuities 
of abrupt changes in the surface coverage of the CNTs. 


We concentrate here on the situations where there is an excess of polymers and maximal cov¬ 
erage of the CNT surface is expected. Binding energies for the maximal coverage of CNTs by 
polymers are shown in Fig. [^(b). Note the discontinuities in the binding energy that are due to 
a sudden change in the number of polymers needed for maximal surface coverage.!^ The posi¬ 
tions of the discontinuities are different from those reported previously in the literature^ 1 5 1 1 8|20|23 1 
because our MD simulations lead to different surface coverages than the geometry optimizations 
performed in those works .1^ These discontinuities have a direct relation to the diameter preference 
of polymers, as will be discussed below. 

Polymer-assisted dispersion as a competition between adsorption and bundling: The bind¬ 
ing energy ^lone cannot explain the selectivity of the polymer adsorption method, 

because its magnitude simply increases with the diameter (see Fig. |^. This would imply that 
large-diameter CNTs are more easily dispersed than small-diameter ones. However, this is in 
clear contrast to the experimental observations discussed above.^^^^^ The binding energy between 
polymer-wrapped CNTs could explain well the polymer-assisted dispersion of CNTs in certain sol¬ 
vents but not the selectivity on CNTs. The key factor for understanding the selection mechanism 
is competition between the bundling of CNTs on the one hand and the adsorption of polymers on 
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the CNT surface on the other (see Fig. [^. This reasoning is based on the observation that CNT 
dispersions in toluene without polymers are not stable and the CNTs eventually rebundle. For this 
competition to take place, the initial state to be considered is a transition state consisting of individ¬ 
ual polymers and individual CNTs. This transition state is experimentally realized by sonication, 
an integral work step of all selection methods. Therefore, the selectivity of CNTs is determined by 
the dijference between the binding energy for CNT bundling and the binding energy for polymer 
adsorption, which reads 

Apbinding _ ^binding _ rebinding 

^ ~ ^CNT-Polymer ^CNT-CNT 



Figure 4: The energetics of CNT bundling: (a) Weighted average Ej (Eq. ??) of the pair binding 
energy of the interaction between a CNT of a given diameter and another CNT, arbitrarily selected 
from the sample, (b) Average number of neighbors Navg of a CNT with a given diameter in a bundle 
with mixed diameters (fractional numbers are a result of the non-uniform diameter distribution), (c) 
Binding energy £’c)vr-^A^r ^ EiNavg of CNT bundling. The variation of £’c)vr-^A?r with the diameter 
follows the same trend as competition between adsorption and 

bundling leads to selectivity. 


Binding energy of CNT bundles: As-produced CNTs are normally a mixture of different 
diameters. This poly disperse nature makes a direct simulation of bundling computationally very 
expensive. To overcome the difficulties, we first calculated the average pair binding energy Ej of 
the interaction between a CNT of a given diameter and another CNT, arbitrarily selected from a 
sample of mixed CNTs. It reads 

E, = '£wiE,j, (3) 

j 

where Ejj is the pair binding energy between CNT species i and j, and wj is the population weight 
(abundance) of CNT species j in the sample. As shown in Fig. |^(a), the magnitude of Ei increases 
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with the CNT diameter, due to the increase in the contact area between CNTs. Next, we estimate 
the average number of neighbors Navg of a CNT in bundles. A simple approach would be to ignore 
the polydispersity and assume that all CNTs have just six neighbors. But our method is to consider 
the surface of a CNT of a given diameter to be covered by CNTs having the average diameter of 
the considered sample, i.e. HiPco or PLV. Therefore, the average number of neighbors can be a 
non-integer. The estimates of Navg for two CNT samples are presented in Fig.|^(b). Finally, the 
binding energy for CNT bundling calculated as 


j^binding _ 

^CNT-CNT ~ 


EiNavg- 


(4) 


As shown in Fig. (c) for both samples, the magnitude of the binding energy of CNT bundles 
increases with the CNT diameter. 

Binding energy difference and diameter selectivity: As discussed already, for both CNT 
bundles and CNT-polymer complexes, the magnitude of the binding energy increases with the tube 
diameter. In the binding energy difference the two trends nearly compensate for each 

other and only their competition leads to the preference for certain diameters, which are reflected 
in the location of the minima of in Fig. 

Consider first the adsorption of PFO on HiPco CNTs in Fig Except for the CNT with the 
smallest diameter, ^binding increases with tube diameter. Since the number of polymers needed 
for the maximal surface coverage of the CNTs changes from three to four, abruptly 

changes at 0.83 nm (see inset of Fig. |^(b)), causing to have a minimum at about the same 

diameter. The behavior of explains (i) the preference of PFO to disperse HiPco CNTs in 

the diameter range 0.8-0.95 nm, a fact that has been repeatedly reported by different groups,l ^b8 | 23 | 
and (ii) why CNTs with diameter smaller than 0.8 nm are not well-dispersed by PFO. These two 
insights explain the dominance of (8,6) CNTs (d = 0.95 nm) in HiPco CNT dispersions and the 
elimination of (6,5) CNTs (d = 0.75 nm) in CoMoCAT CNT dispersions. 

For PFH, with side-chains two carbon atoms shorter than PFO, more polymer chains are needed 
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to cover the surfaee of a CNT. Therefore, the diseontinuity in /^binding upshifted to 1.03 nm. 
This explains why, for HiPeo CNTs using PFH instead of PFO, the dominant CNTs in the disper¬ 
sion beeome (8,7) (d = 1.02 nm) and (9,7) (d = 1.09 nm) (see Fig. lb of Nish et alJ^. 

For the eopolymer PFD-Py, the minimum of is at about 1.25 nm. This agrees with re- 

eent experimental findings that, for HiPeo CNTs, PFD-Py prefers to disperse CNTs with diameters 
of about 1.23 nm (see Fig. In of Berton et al^). 

The of PFH-A inereases eontinuously in the eonsidered diameter range (0.8-1.4 nm) 

and no minimum is diseemible. Mistry et al. performed a systematie study on the seleetivity of 
PFH-A on CNTs synthesized via laser vaporization of graphite at different temperatures and found 
that it always prefers to disperse CNTs with the smallest diameters in the sample.^^^ The absenee 
of a minimum in that range is again eonsistent with the experiments even though the simulations 
are based on HiPeo CNTs. Further diseussions on the seleetivity of PFH-A ean be found in the 
Supporting Information. 



Figure 5: Diameter seleetivity as eompetition between CNT bundling and polymer adsorption. The 
diameter preferenees of speeifie polymers for HiPeo CNTs in our simulations are defined by the 
minima of the eorresponding binding energy differenee (Eq. ??). Th ey are in exeellent 

agreement with experimental results indieated by the shadowed regions 


To summarize, for the polymers PFO, PFH and PFD-Py, the minima of the binding energy 
differenee mateh perfeetly the experimentally reported diameters that are dominantly 

dispersed by those polymers. This exeellent agreement strongly suggests that the meehanism of 
diameter seleetivity is a eompetition between CNT bundling and polymer adsorption. 

It is interesting to note that the sign of is negative for PFO, PFH and PFD-Py. This 
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indicates a preference for the formation of CNT-polymer adsorption complexes over CNT-CNT 
bundling. Therefore, for long sonication times, SWNTs of all diameters can be dispersed in prin¬ 
ciple. With increasing sonication time, the amount of dispersed CNTs will increase and the se¬ 
lectivity will gradually diminish. Therefore, an optimal sonication time should be experimentally 
determined, providing a compromise between yield and purity. Positive values of for 

PFH-A mean that the rebundling should happen more frequently than adsorption, which implies a 
potential lower yield of the dispersion process using PFH-A. 

Note also that, due to the possibility of partial adsorption and other "imperfect" packing con¬ 
figurations the transition in Fig. 3(b) may turn out to be not so abrupt, see also the radial distribu¬ 
tion functions shown in Fig.S6 and the corresponding discussions in the Supporting Information. 
Therefore, in experiments a range of diameters is often selected by a certain polymer. 

One can also view the sonication-assisted dispersion process as a reversible reaction, 

CNT@CNT ^PFO^PFO@CNT. (5) 

In this language, the initial configuration of isolated CNTs and polymers corresponds to a transition 
state, which is achieved with the aid of ultrasonication treatment adsorption of polymers 
on the CNT surface and the bundling of CNTs are the rate-determining steps for forward and 
backward reactions, respectively. The two binding energies are the (negative) activation energies 
for reactions in the two directions. For this reversible reaction, the binding energy difference only 
estimates energetic contributions to the reaction rates, neglecting entropic contributions, reaction 
orders and the concentrations of the reactants. 

The focus of the current study is on the diameter selectivity of CNTs by aromatic polymers. For 
the chirality selection, the match/mismatch between the atomic structures of polymers and CNTs 
will be quite crucial. The popular implementations of van der Waals’ interaction, as used here, 
were found unsuitable for the purpose and the anisotropic intermolecular potentials turn out to be a 
better alternative.^ For the sorting of CNTs with respect to electronic properties, ab initio quantum 
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simulations with the electronic interactions being included would be more appropriate. All these 
issues deserve their own separate publications. The mentioned success of our simple energetic 
model implies that the key factors determining the diameter selectivity of (semi-conducting) CNTs 
were properly represented. The model can certainly be further improved by including: (i) entropy 
factors for a proper estimate of Gibbs free energy, which is of direct relevance to the reaction 
kinetics, and (ii) the effect of explicit solvent. Calculations with explicit solvents and estimates of 
entropic contributions are provided in the Supporting Information. 

In conclusion, we explain the diameter selectivity of polymer adsorption methods to be the 
result of a competition between the bundling of CNTs and the adsorption of polymers on the CNT 
surface. The preference of certain diameters corresponds to local minima of the binding energy 
difference between these two processes. Such minima occur due to abrupt changes in the CNT’s 
coverage with polymers at certain diameters. For all tested polymers including two homopolymers 
of polyfluorene with different side-chain lengths and two copolymers with anthracene or pyridine 
groups, our simulation results are in excellent agreement with the experimental findings regrading 
the diameter selectivity. Interestingly, even the influence of a fine-tuning of side-chain length on 
the selectivity was correctly captured in our method. Our insights resolve a long-standing con¬ 
troversy regarding the understanding of CNT selection schemes and are important for the further 
development of dispersion/extraction methods, i.e., they enable MD simulations to be used for the 
screening of polymer candidates, tailoring of polymer structures, and obtaining further scientific 
insights. The proposed mechanism is general enough to be valid for other (sonication-aided) dis¬ 
persion processes, for instance, the exfoliation of layered materials,!^ and the dispersion of CNTs 
by DNAs and mononucleotides E5|47 | 49|50| 

Methods 

The adsorption of polymers on single-walled CNTs and the bundling of CNTs were studied with 
classical molecular dynamics (MD) simulations by using the CP2 kI^ and Gromacs packages.® 
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MD simulations were performed in NVT ensemble at T = 300^ using the Nose-Hoover or Langevin 
thermostats. The standard CHARMM foree field parameters for the intra-moleeular interaetion^ 
were benehmarked against the density funetional method (DFT) MD simulations with Grimme 
dispersion eorreetions DFT-D3®^ and the BLYP exehange-eorrelation funetional^ in CP2K and 
elassieal MD simulations with the MM3 foree field® using the Tinker paekage.® The torsion an¬ 
gle parameter, deseribing amongst others the twist of the baekbones of polymers, was modified 
to mateh the results of DFT first prineiples MD simulation. The inter-moleeular interaetions in- 
elude an eleetrostatie part due to partial eharges on atoms and a dispersion foree part modeled by 
a Lennard-Jones potential as usual in standard CHARMM foree field implementations.® 

For the adsorption of polymers on the CNT surfaee, the polymer baekbones were initially 
aligned parallel to the CNT axis. Multiple ehains of polymers were arranged in an n-fold symmetrie 
strueture surrounding the tube. The initial distanee between the baekbone and the CNT surfaee 
was set to 1 to 1.5 nm depending on the CNT diameter and the number of polymers. For the 
binding energy ealeulation of CNT pairs, the two CNTs were plaeed in parallel with the initial 
distanee between the surfaees of 0.6 nm. The time step for the integration of Newton’s equation of 
motion was 1 fs. The duration of MD simulations ranged from 1 ns to 20 ns. For the ealeulation of 
thermodynamie averages, the equilibration time, ranging between 0.2 and 2 ns, was not eonsidered. 
To eheek the stability of self-eonstrueted helieal adsorption struetures, geometry optimizations 
were performed using the CP2K paekage. The eriteria of eonvergenee are 3 x 10^^ Bohr for the 
geometry ehange and 4.5 x 10^^ Hartree/Bohr for the ehange in the foree. 

It is known that the solvent plays an important role in the seleetive dispersion of CNTs by 
polymers.® However, here we are interested in studying the effeet of different polymers in eom- 
bination with the same weakly polarized solvent, toluene. Moreover, our tests showed that the 
explieit inelusion of solvent moleeules of toluene in MD simulations did not ehange the struetures 
of adsorption but signifieantly slowed down the dynamies of the adsorption proeess. Therefore, to 
enable MD simulations within reasonable times, solvent moleeules were not explieitly ineluded in 
our produetion runs. Our additional MD simulations with explieit solvents showed that, the ad- 
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sorption configurations of the polymer baekbones and sidechains, whieh are in contact with CNT 
surface, hardly change with the inclusion of solvents. Only the sidechains, which are not in eon- 
taet with CNT surfaee, tend to point outwards to the solvent instead of folding baek and aligning 
along CNT surfaee as in vaeuum. The seeond-layer of polymers moves a bit away from the CNT 
surfaee. Therefore, with the inelusion of solvents the values of binding energies may ehange by 
some degree but the surfaee eoverage of CNTs and the related positions of the abrupt changes in 
Fig.3(b) will be unaffected. Further details on the effeet of explieit solvent ean be found in the 
Supporting Information. Studies showed that a PFO octamer already has the same selectivity as 
a PFO polymer. Furthermore, the stability of the oligomer-CNT eomplex inereases strongly with 
the ehain length of the oligomer .1^ To meet the eapaeity of the available eomputing resourees in 
our simulation, we used 30-nm-long polymer ehains, whieh eonsist of 32, 22 and 30 monomers 
of PFO/PFH, PFD-Py and PFH-A, respeetively. CNT segments of length between 30 and 36 nm 
were used, their length varying with the ehirality. 

To obtain the binding energy of a polymer-CNT eomplex, three MD simulations were per¬ 
formed for the adsorption eomplex, the isolated CNT and the polymer, respeetively. The mean 
value of the potential energy was determined from the eorresponding trajeetories and the binding 
energy was then ealculated. This procedure is different from most oases in the literature where the 
binding energy was oaloulated from {T = 0 K) single-point calculations of the optimized structures. 

It is worth pointing out that, for the adsorption of polymers on a CNT surfaee, the binding 
energy ean be measured per unit length of a polymer ehain, or per unit length of CNT oovered 
oompletely by polymers. The former desoribes how hard it is to remove a polymer ehain from the 
CNT surface while the latter is suitable for characterizing the competition for the adsorption on a 
CNT surfaee. 
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